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ABSTRACT 

We  have  considered  stationary  electric  and  magnetic  fields  of  fractal  objects  with  the 
help  of  differintegration  methods.  The  equations  of  fractal  electrostatic  and 
magnetostatic  potential  at  integro-differential  form  are  proved  and  introduced.  The 
offered  ideas  can  be  useful  for  modification  numeric  fields  modelling  techniques  to 
solve  electromagnetic  problems  for  real  objects  with  material  structure  irregularities 
consideration. 


INTRODUCTION 

There  are  a number  of  modern  numeric  electromagnetic  fields  modelling  techniques. 
Finite-Difference  Time-Domain  (FDTD),  Method  of  Moments  (MoM),  and  Finite 
Element  Method  (FEM)  are  the  most  popular  and  useful  methods.  The  general  essence 
of  these  methods  consists  in  division  of  the  approached  ideal  model  of  real  physical 
object  into  elementary  components  the  fields  for  which  can  be  found  on  the  basis  of  the 
classic  Maxwell’s  theory.  But  the  classic  electromagnetic  theory  in  which  the  object’s 
geometry  definition  is  based  on  concepts  of  point,  line  and  plane  (ideal  objects)  become 
impenetrable  when  the  explanation  of  the  field  distribution  of  roughness  fractal  surfaces 
is  needed.  The  examples  of  differintegration  methods  application  for  such  complex 
electromagnetic  problems  were  originally  shown  at  [1]. 

MAIN  PART 


Based  on  idea,  that  the  fractal  sets  are  adequate  geometrical  model  for  irregular  contours 
and  surfaces.  Let  us  consider  prefractal  covering  compact  set,  which  is  our  proposed 
model  of  a contour  with  a current,  as  a limit  of  monotonously  growing  sequence  of 
covering  compact  sets  with  the  corresponding  sequence  of  rising  contours 
/ c /,  c ...  cz  c ...  /„ . The  current  near  by  /„  contour  can  be  found  as: 

J x /„_!  1 1 


^X>  = - Adln~\  \dln_ 2 \p(x')dx'  . 


A{l,n)  * 

The  obtained  repeated  integral  may  be  summarize  by  Riemann's-Liouville's  definition 
of  fractional  integral 


LDaxp){x)  = 


1 


p(x') 


T(a)  l (x 


■x'ta 


dx' . 


(1) 
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As  a result  of  usage  fractal  contour  (surface)  covering  by  compact  sets,  we  have  reduced 
a problem  to  construction  a smoothing  Hausdorfs  measure  on  a physically  prefractal 
layer  with  differintegration  of  an  uniform  electric  charge  density  on  a fractal  point  set 
projection  to  a smooth-faced  segment  [2]. 

Using  of  a -order  differintegral  definition  of  fractal  electric  charge  density  pa  ( r') 
located  at  r'  Eq.  (1)  enables  to  formulate  the  equations  of  classical  electrostatics  in 
terms  of  a -characteristics.  In  analogy  with  the  fractal  electrostatic  case,  we  may  define 
fractal  electric  current  density  jft(r)  for  magnetostatic  field  analysis.  So  the  laws  of 
electrostatics  and  magnetostatics  for  fractal  objects  can  be  summarized  in  two  pairs  of 
time-independent,  uncoupled  vector  differintegral  equations,  namely  the  equations  of 
fractal  electrostatics 

V • E"(r)  = _ (2) 

V • E"  (r)  = 0 , 

and  the  equations  of  fractal  magnetostatics 


VxH"(r)  = //0j“(r).  (3) 

V x H " ( r)  = 0 . 

The  electrostatic  field  E"(r)  is  irrotational  and  it  may  be  expressed  in  terms  of  the 
gradient  of  a scalar  field.  If  we  denote  this  scalar  field  by  -^"(r),  we  get 
Eft(r)  = -V^"(r) . Taking  the  divergence  of  this  and  using  Eq.  (2),  we  obtain  Poisson’s 
equation 

p"(r) 


V~(j)a  (r)  = -V  • E"  (r)  = -- 


The  solution  of  Eq.  (4) 


Fir)  = 


Pa(  r) 


4 US, 


civ'  + c , 


(4) 


(5) 


o r 


where  the  integration  is  taken  over  all  source  points  r'  at  which  the  charge  density 
pa ( r')  is  non-zero  and  c is  an  arbitrary  quantity  which  has  a vanishing  gradient.  The 
scalar  function  </>n(r)  in  Eq.  (5)  above  is  called  the  fractal  electrostatic  scalar  potential. 
Consider  the  equations  of  magnetostatics  Eq.  (3)  we  got  fractal  magnetostatic  vector 
potential  with  definition  from  EC (r)  = VxA“ (r)  as 


(6) 


A"(r)=fLIr^'/v'+a(r)- 

An  * I r - r ( 

where  a(r)  is  an  arbitrary  vector  field  whose  curl  vanishes. 

CALCULUS  EXAMPLES 

As  an  example  of  using  obtained  expressions  (Eq.  4-6)  at  FEM  methods  applications  let 
us  consider  following  problem 


V‘0  = - 1 inside  B,  Da <f>  =0, 


(7) 
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where  A border  is  a circle  defined  as  {x  ■ cos(l);  y ■ sin (/);  0 < t < 2n\  and  B circle  is 
infinity  approximation  defined  as  {x  • 5cos(/);  y • 5sin(/);  0 < / < 2k). 


= C/(x-b)'~a,  \ <a<2. 


Fractal  boundary  condition  Da<f> j =C  translates  into  (f>{a) 

The  problem  (7)  was  solved  for  b = - 5,  C = 1 parameters  by  finite  element  method 
with  the  help  of  FreeFEM+  software.  The  obtained  potential  surface  for  a - l .7  is 
shown  in  Fig.  1.  We  demonstrate  the  XOZ  -plane  cut  of  potential  distribution  for  a - 
parameter  varied  from  1 up  to  2 in  Fig.  2._ 


Fig.  1.  3D  potential  distribution  for 

a = 1.7 


Fig.  2.  Potential  distribution  for  various  a - 
parameter 


The  problem  (7)  with  a = 1 and  a - 2 parameters  of  fractal  boundary  condition  is 
equivalent  of  Dirichlet. 

More  interesting  that  at  a = 1.5  calculation  result  is  equal  to  solving  of  mixed  problem 
with  Dirichlet  and  Neumann  boundary  conditions  like  (dtp/d n)|^  +<t\4  = x . 

CONCLUSIONS 


The  application  of  fractional  calculus  enables  to  formulate  the  electromagnetic  potential 
equations  in  terms  of  a -characteristics.  Formally  proposed  equations  agree  with  the 
classical,  at  imposing  on  a -characteristics  additional  boundary  conditions. 

Numeric  decision  of  typical  problem  of  potential  theory  has  shown  that  definition  of 
boundary  conditions  in  the  integro-differetial  form  generalizes  classical  boundary 
conditions  and  is  adequate  for  fractal  surfaces  with  various  levels  of  roughness.  The 
offered  ideas  assumes  to  remove  difficulties  in  investigations  of  singular  distributions 
and  can  be  useful  for  modification  numeric  fields  modeling  techniques  to  solve 
electromagnetic  problems  for  real  objects  with  structure  irregularities  consideration. 
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